Nonlinear turbulent magnetic diffusion and effective drift velocity of large-scale 
magnetic field in a two-dimensional magnetohydrodynamic turbulence 
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We study a nonlinear quenching of turbulent magnetic diffusion and effective drift velocity of large- 
scale magnetic field in a developed two-dimensional MHD turbulence at large magnetic Reynolds 
numbers. We show that transport of the mean-square magnetic potential strongly changes quenching 
of turbulent magnetic diffusion. In particularly, the catastrophic quenching of turbulent magnetic 
diffusion does not occur for the large-scale magnetic fields B 3> B cq /vRm when a divergence of the 
flux of the mean-square magnetic potential is not zero, where B oq is the equipartition mean magnetic 
field determined by the turbulent kinetic energy and Rm is the magnetic Reynolds number. In this 
case the quenching of turbulent magnetic diffusion is independent of magnetic Reynolds number. 
The situation is similar to three-dimensional MHD turbulence at large magnetic Reynolds numbers 
whereby the catastrophic quenching of the a effect does not occur when a divergence of the flux of 
the small-scale magnetic helicity is not zero. 

PACS numbers: 47.65.Md 



I. INTRODUCTION 



The magnetic fields of the Sun, solar type stars, galax- 
ies and planets are believed to be generated by a dynamo 
process due to the simultaneous action of the a effect (the 
helical motions of turbulence) and differential rotation 
(see, e.g., [H, 0, H, 0, [E @|)- The kinematic stage of the 
mean-field dynamo, i.e. the growth of a weak mean mag- 
netic field with negligible effect on the turbulent flows, 
is well understood, while the nonlinear stage of dynamo 
evolution is a topic of intensive discussions (for reviews, 
see 0, B El). The most contentious issue is the ques- 
tion of the equilibrium magnetic field strength at which 
dynamo action saturates. In particular, the problem of 
catastrophic quenching of the a effect in a developed 
three-dimensional magnetohydrodynamic (MHD) turbu- 
lence with large magnetic Reynolds numbers has been 
intensively discussed in astrophysics and magnetoh ydr o- 
dynamics during last years (see, e.g., [13, El, [H, [HI])- 
The catastrophic quenching implies very strong reduc- 
tion of the a effect during the growth of the mean mag- 
netic field so that the dynamo generated magnetic field 
should be saturated at a very low level. However, this is 
in contradiction with observations of the magnetic fields 
of the Sun, stars and galaxies. 

In a two-dimensional MHD turbulence with imposed 
large-scale magnetic field at large magnetic Reynolds 
numbers, the catastrophic quenchingcan occur for tur- 
bulent magnetic diffusion (see, e.g., [1(1 [3 )• I n particu- 
lar, small-scale magnetic fluctuations strongly affect the 
large-scale magnetic field dynamics even for very weak 
mean fields. This causes a strong reduction of turbu- 



lent magnetic diffusion 1 101 ] . This conclusion is based on 
Zeldovich theorem [l6| . In a two-dimensional MHD tur- 
bulence energy is transferred from large-scale stirring to 
small scales and dissipated due to an Alfvenized cascade, 
whereby eddy energy is converted to Alfven wave energy 
(see, e.g., [13, [Hj]). The above discussed quenching is 
caused by the tendency of the mean magnetic field to 
Alfvenize the turbulence. 

A principal difference between two-dimensional and 
three-dimensional MHD turbulence is related to differ- 
ent integral of motions for these kind of turbulence. In 
particular, square of total (small-scale and large-scale) 
magnetic potential is conserved in two-dimensional MHD 
turbulence, while total (small-scale and large-scale) mag- 
netic helicity is conserved in three-dimensional MHD tur- 
bulence. The magnetic helicity and the a effect can be 
positive and negative, while the squared magnetic po- 
tential is only positive. A comprehensive comparison be- 
tween two-dimensional and three-dimensional MHD tur- 
bulence has been performed in [3, [l5j . 

It has been recently recognized [l9|, US] that in three- 
dimensional MHD turbulence the catastrophic quench- 
ing of the a effect does not arises when a divergence 
of the flux of magnetic helicity is not zero (see also 
@)IHl[13l)- I n the present study we show that in a devel- 
oped two-dimensional MHD turbulence with large mag- 
netic Reynolds numbers Rm, a non-zero divergence of 
the flux of the mean-square magnetic potential strongly 
changes a balance in the equation for these fluctuations 
and results in that the catastrophic quenching of turbu- 
lent magnetic diffusion does not occur for the magnetic 
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where B cq is the equipartition 
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mean magnetic field determined by the turbulent kinetic 
energy. 

This paper is organized as follows. In Sec. II we formu- 
late the governing equations, the assumptions, the pro- 
cedure of the derivations. In Sec. Ill we determine the 
nonlinear turbulent magnetic diffusion coefficients and 



2 



the nonlinear drift velocities of the mean magnetic field 
in a developed two-dimensional MHD turbulence. Fi- 
nally, we draw conclusions in Sec. IV. In Appendix A we 
perform the derivation of the nonlinear turbulent mag- 
netic diffusion and the nonlinear drift velocities of the 
mean magnetic field and in Appendix B we present the 
nonlinear functions used in Sec. Ill and their asymptotic 
formulas. 



II. GOVERNING EQUATIONS AND THE 
PROCEDURE OF DERIVATION 

Let us consider a developed two-dimensional MHD tur- 
bulence with large hydrodynamic and magnetic Reynolds 
numbers. We study nonlinear quenching of the turbu- 
lent magnetic diffusion and the effective drift velocity 
of the magnetic field. We use a mean field approach 
whereby the velocity, pressure and magnetic field are sep- 
arated into the mean and fluctuating parts. In a two- 
dimensional MHD turbulence the mean magnetic field is 
B = Vx [A(x, y) e], where A{x 1 y) is the mean magnetic 
potential and e is the unit vector perpendicular to the 
plane of the two-dimensional MHD turbulence, i.e., it is 
directed along z-axis. The equation for the evolution of 
the mean magnetic potential for an incompressible tur- 
bulent flow with a zero mean velocity reads: 



dA 
~dt 



+ div (u a) = tjAA 



(1) 



where u are the velocity fluctuations and r\ is the mag- 
netic diffusion caused by an electrical conductivity of a 
fluid. The mean electromotive force is £ z = (u x b) 2 = 
— divT a , where the spatial flux of magnetic potential 
r a = (us) and magnetic fluctuations b = Vx[a(x,y)e] 
are described by the fluctuations of the magnetic poten- 
tial a{x,y). The mean electromotive force f^(B) in a 
two-dimensional MHD turbulence is given by: 

£ Z (B) = {[V cff (B)xB] J -7 7u (B)(VxB),}e l , (2) 

where the nonlinear turbulent magnetic diffusion r]ij (B) 
and the nonlinear effective drift velocity V cff (B) of the 
mean magnetic field are determined in Sec. III. 

In order to derive equations for the nonlinear turbu- 
lent magnetic diffusion and the nonlinear effective drift 
velocity of the mean magnetic field in a two-dimensional 
MHD turbulence we use a procedure outlined below (see 
Appendix A for details). This procedure is similar to 
that used in [23[ for a study of a three-dimensional MHD 
turbulence. We use equations for fluctuations of velocity 
and magnetic field 



du(t, x) 



db(t,x) 



Vp 

P 



I 

1 

fip- 

F. 



[(b • V)B + (B ■ V)b] 



(B • V)u — (u • V)B + b 



N 



(3) 
(4) 



where p is the fluid density, F is a random external stir- 
ring force, u N and h N are the nonlinear terms which 
include the molecular dissipative terms, p are the fluc- 
tuations of total (hydrodynamic and magnetic) pressure. 
Hereafter we omit the magnetic permeability of the fluid 
p, and include p^ 1 / 2 in the definition of magnetic field, 
we also omit the density p of incompressible fluid and 
include p 1 / 2 in the definition of velocity field. We rewrite 
Eqs. (J3]) and ((4]) in a Fourier space and derive equations 
for the two-point second-order correlation functions of 
the velocity fluctuations (uiUj), the magnetic fluctua- 
tions (bibj) and the cross-helicity (biUj). The equations 
for these correlation functions are given by Eqs. (|A1[) - 
(|A3]) in Appendix A. 

The second-moment equations include the first-order 
spatial differential operators Af applied to the third-order 
moments M^ ITI \ A problem arises how to close the sys- 
tem, i.e., how to express the set of the third-order terms 
MM^ 111 ^ through the lower moments M^ 11 ' (see, e.g., 
24, 25, 26]). We use the spectral r approximation which 
postulates that the deviations of the third- moment terms, 
AfM^- 111 ^ (k), from the contributions to these terms af- 
forded by the background turbulence, NM^ III,0 '(k), are 
expressed through the similar deviations of the second 
moments, M( 77 )(k) - Af( 7/ '°)(k): 



A/"M (/7/) (k) 



Af M (///,o) (k) 

-J-[M(")(k) 
T{k) 



M^°)(k)] , (5) 



(see, e.g., [H H [23, [H HH), where r(k) is the 
scale-dependent relaxation time, which can be identi- 
fied with the correlation time of the turbulent velocity 
field, and the quantities with the superscript (0) corre- 
spond to the background turbulence. A justification of 
the r approximation for different situations has been per- 
formed in numerical simulations and analytical studies in 

[1 HIM 11 III, [11. 

Next, we split all second-order correlation functions, 
M^ 11 *, into symmetric and antisymmetric parts with re- 
spect to the wave vector k. We assume that the char- 
acteristic time of variation of the mean magnetic field 
B is substantially larger than the correlation time r(fc) 
for all turbulence scales. This allows us to get a sta- 
tionary solution for the equations for the second-order 
moments, M( n \ Wc use a model of the background 
anisotropic and inhomogeneous two-dimensional MHD 
turbulence determined by Eqs. (|A16|) - (|A17|) in Appendix 
A. 

In this study we consider an intermediate nonlinearity 
which implies that the mean magnetic field is not enough 
strong in order to affect the correlation time of turbulent 
velocity field. The theory for a very strong mean mag- 
netic field can be modified after taking into account a 
dependence of the correlation time of the turbulent ve- 
locity field on the mean magnetic field. 

Using the solution of the derived second-moment equa- 
tions, wc determine the mean electromotive force, E, — 
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£imn /(frn w m )k dk (see Appendix A for details), where 
Sijk is the fully antisymmetric Levi-Civita tensor. This 
procedure allows us to determine the nonlinear turbulent 
magnetic diffusion and the nonlinear effective drift veloc- 
ity of the mean magnetic field in a two-dimensional MHD 
turbulence. 



III. TURBULENT TRANSPORT 
COEFFICIENTS 

The derivation outlined in Sec. II yields the nonlin- 
ear turbulent magnetic diffusion of the mean magnetic 
field. In particular, in order to determine the nonlinear 
turbulent magnetic diffusion rjij (B) we use an identity: 
Vij = (£ik P bjkp +£ jkpbjk p)/4:, where the tensor b ljk is de- 
termined by Eq. (|Al9j) in Appendix A. The nonlinear 
turbulent magnetic diffusion coefficient along the mean 
magnetic field, r\ B , and the cross-field turbulent magnetic 
diffusion coefficient, 77^, are given by: 



V± = T [(u 2 )( 0) -<b 2 ) (0) W), 



(6) 
(7) 



where t = Io/uq and w = y/ (u 2 )(°) is the characteris- 
tic turbulent velocity in the maximum scale of turbulent 
motions ^o- The quantities with the superscript (0) cor- 
respond to the background turbulence. The functions 
^(/3), x f'i(/3) and their asymptotic formulas are given in 
Appendix B, f3 = 4 (B/B cq ) and B cq = ^(u 2 )(°) is the 
equipartition field. More general equations for r\ B and 
rj ± in the case of an anisotropic background turbulence 
are given by Eqs. (|A"20)) and (|A"2i"|) in Appendix A. It 
follows from Eqs. © and ([7]) that in the case of Alfvenic 
equipartition, (u 2 )(°) = (b 2 )(°), the nonlinear turbulent 
magnetic diffusion vanishes. 

The nonlinear turbulent magnetic diffusion depends on 
a flux of mean-square magnetic potential. This flux can 
change properties of the quenching of the cross-field tur- 
bulent magnetic diffusion. Indeed, let us determine the 
parameter e = (b 2 )(°)/(u 2 )(°) using budget equation for 
the evolution of the mean-square magnetic potential (a 2 ): 

S(a 2x 



divF A = 2r^B 2 - 2?7(b 2 



(8) 



(ua 2 ) — 77V (a 2 ) determines the 



dt 

where the flux F'' 
transport of (a 2 ). The first term 2?7 i B 2 in the right 
hand side of Eq. ([5]) describes a production of the mean- 
square magnetic potential (a 2 ), while the term — 2?7(b 2 ) 
determines the resistive dissipation of (a 2 ). In the ab- 
sence of the flux of the mean-square magnetic potential, 
F A — 0, Eq. implies the catastrophic quenching of 
the cross-field turbulent magnetic diffusion. In partic- 
ular, in a steady-state Eq. ([5]) reads rj ± = rj(b 2 )/B 2 . 
Since the magnetic energy is less than the kinetic energy, 



< (u 



2\(0) 



we get 

!k< 1 

Vt Rm(B/B cq )2 



(9) 



where r\ T — Iqu /2 and Rm = u l /i] is the mag- 
netic Reynolds number. This estimate implies a strong 
quenching of the cross-field turbulent magnetic diffusion 
with increasing Rm (B/B cq ) 2 due to Alfvenization of tur- 
bulence by tangling of a weak mean magnetic field by 
velocity fluctuations flij ]. 

Situation is drastically changed when div F A ^ 0. In- 
deed, Eq. ([5]) is not closed because it depends on the 
magnetic energy (b 2 ). The energy of magnetic fluctu- 
ations (b 2 ) can be determined in the same way as we 
derived the cross-helicity tensor. In particular, (b 2 ) is 
obtained from Eq. (|A8[) given in Appendix A, after the 
integration in k-space. The result is given by 



, 2 \(°)h - 



[i-<K/?)] + (b 2 )(°)[i + </>(/?)] 



(10) 



where the function </>(/?) and its asymptotic formulas are 
given in Appendix B. More general equation for (b 2 ) for 
anisotropic background turbulence is given by Eq. (|A22[) 
in Appendix A. 

Equation ([5]) allows us to determine the energy of 
magnetic fluctuations of the background turbulence self- 
consistently. In particular, combining Eq. (|10p with the 
steady-state solution of Eq. ((8|) we determine the pa- 
rameter e = (b 2 )(°V(u 2 ) W [see, e.g., Eq. (jX23| in Ap- 
pendix A for anisotropic background turbulence] . When 
B 3> B cq / v / Rm, the parameter e is given by 



e = 1 



divF A 

4t7 t S 2 *(/3) ' 



(11) 



Therefore, Eqs. (|6|), (O and (fTTjl yield the nonlinear tur- 
bulent magnetic diffusion in two directions: 



div F J 



r h 



2B 2 V 
divF A 



2B 2 



(12) 
(13) 



where r) B is the nonlinear turbulent magnetic diffusion 
along the mean magnetic field and is the cross-field 
nonlinear turbulent magnetic diffusion. Remarkably, 
Eq. (fT3)l can be obtained directly from Eq. ((8|) written 
in a steady-state if we neglect the resistive dissipation 
term — 2?7(b 2 ) in the right hand side of Eq. ©. 

In order to determine the parameter e we use the 
steady-state solution of Eq. (JSJ). However, the steady- 
state solution of this equation exists not for all values of 
the mean magnetic field. Indeed, let us plot in Fig. 1 the 
function ^ifi) for the exponent of the energy spectrum 
of the background turbulence q = 5/3. At B — > 0.18 B cq 
the function \&(/3) tends to zero (see Fig. 1). In the range 
B > 0.18 B eq the steady-state solution of Eq. (jHJ) does not 
exist. The turbulent magnetic diffusion should be pos- 
itive, which implies that divF" 4 > 0. Therefore, when 
divF A < there is no steady-state solution of Eq. ([5]) 
for B > 0.18 £> C q as well. More detailed discussion of this 
facet is given in Appendix A after Eq. (j A26|) . 
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0.2 0.4 0.6 0.8B/B 
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FIG. 1: The function *(B/B cq ) for q = 5/3. 



In inhomogeneous turbulence there are also turbulent 
diamagnetic and paramagnetic effects. In particular, 
an inhomogeneity of the velocity fluctuations leads to a 
transport of mean magnetic flux from regions with high 
intensity of the velocity fluctuations (turbulent diamag- 
netism, see, e.g., [HEH)- On the other hand, an inhomo- 
geneity of magnetic fluctuations due to the small-scale 
dynamo causes turbulent paramagnetic velocity, i.e., the 
magnetic flux is pushed into regions with high intensity of 
the magnetic fluctuations (see, e.g., [35ll36l|). In order to 
determine the nonlinear turbulent diamagnetic and para- 
magnetic drift velocities V cff (B) of the mean magnetic 

field, we use an identity: vjf^ = ej^uaij/2, where the 
tensor is determined by Eq. (|A18|1 in Appendix A. 
The inhomogeneities of the velocity and magnetic fluctu- 
ations of the background turbulence are characterized by 
the following parameters = V;(u 2 )(°)/(u 2 )(°) and 

A| 6) = V 4 (b 2 )(°V(b 2 )(°). The nonlinear effective drift 
velocity of the mean magnetic field is given by 

V cff = -2?? T [A M - eA( m )]*!(/3) , 



(14) 



where the function M/i and its asymptotic formulas are 
given in Appendix B. When B > B cq /VRm, Eqs. (fTTj) 
and HU) yield 



V 



off 



Hiv F^ 



(15) 

The first term cx AS"' in Eq. (fl~5|) determines the tur- 
bulent diamagnetic drift velocity while the second term 
cx A^™) describes the turbulent paramagnetic drift veloc- 
ity. The last term oc divF' 4 in Eq. (fT5)) determines the 
turbulent diamagnetic drift velocity caused by magnetic 
fluctuations for B < 0.18 B cq . More general equation for 
V for anisotropic background turbulence is given by 
Eq. (|A28|) in Appendix A. 



IV. CONCLUSIONS 

In the present study we investigate nonlinear quench- 
ing of the turbulent magnetic diffusion and the effective 
drift velocity of the magnetic field in a developed two- 
dimensional MHD turbulence at large magnetic Reynolds 
numbers. We elucidate an important role of transport 
of the mean-square magnetic potential which strongly 
changes quenching properties of turbulent magnetic dif- 
fusion. In particular, we show that the catastrophic 
quenching of turbulent magnetic diffusion does not arises 
for the magnetic fields B ^s> B cq /v / Rm for a non-zero 
divergence of the flux of the mean-square magnetic po- 
tential. In this case the quenching of turbulent magnetic 
diffusion is independent of magnetic Reynolds number. 
This is similar to a three-dimensional MHD turbulence 
at large magnetic Reynolds numbers whereby the catas- 
trophic quenching of the a effect does not occur when 
a divergence of the flux of the small-scale magnetic he- 
licity is not zero. Note that in a two-dimensional MHD 
turbulence, the magnetic field may only decay, while in 
three-dimensional MHD turbulence magnetic field may 
grow by dynamo mechanism. 

Note that a quenching of turbulent magnetic diffusiv- 
ity in a 'wavy' magnetohydrodynamic turbulence in two 
dimensions was recently studied in [3^ | . They found that 
the turbulent magnetic diffusivity in the fourth-order 
does not vanish when the magnetic Reynolds number 
tends to infinity. In particularly, the second-order (quasi- 
linear) contribution to the spatial flux of the mean mag- 
netic potential is quenched as Rm _1 , while the fourth- 
order contribution to the flux is independent of Rm. 
This implies that the turbulent magnetic diffusivity is 
not quenched catastrophically in the presence of disper- 
sive waves which can transfer the mean-square magnetic 
potential. These findings are in an agrement with our 
results. 
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APPENDIX A: DERIVATIONS OF THE 
NONLINEAR TURBULENT TRANSPORT 
COEFFICIENTS 

We use equations for fluctuations of velocity and mag- 
netic field written in a Fourier space and derive equations 
for the second moments in two-dimensional MHD turbu- 
lence using a procedure which is similar to that used in 
[23l | for a study of a three-dimensional MHD turbulence. 
In order to exclude the pressure term from the equation 
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of motion (|3|) we determine Vx(Vxu). We also apply 
the two-scale approach, e.g., we use large-scale R, K and 
small-scale r, k variables (see, e.g., [38|). We assume that 
there exists a separation of scales, i.e., the maximum scale 
of turbulent motions is much smaller then the charac- 
teristic scale Lb oi inhomogeneities of the mean magnetic 
field. We derive equations for the following correlation 
functions: /y(k,R) = L(ui;uj), /iy(k, R) = L(bi;bj) 
and <7jj(k, R) = L(bf,Uj), where 

L(a; c)= J (o(k + K/2)c(-k + K/2)) 
x exp (iK-R) dK . 

The equations for these correlation functions are given 
by 

dfM = i^B)^+4 + F ij+ Ar fij , (Al) 



Of 



^ = -i(k.B)*W+Z*+^ 



(A2) 



%j( k ) 
dt 



= i(k-B)[/ i3 -(k) - %(k)] + /* + , 



(A3) 



where hereafter we omit arguments t and R in the cor- 
relation functions and neglect terms ~ 0(V 2 ). Here 
= fl«(k) - fti (-k), F y (k) = (F(k) Wj (-k)) + 

(^(kJF^-k)}, and F(k) = kx (kxF(k))/fc 2 p. The stir- 
ring force F(k) is an external parameter, that determines 
the background turbulence. The source terms ifj , 
and /?■ which contain the large-scale spatial derivatives 
of the mean magnetic field and turbulence are given by 



We use the spectral r approximation which postu- 
lates that the deviations of the third-moment terms, 
AfM^ 11 ^ (k), from the contributions to these terms af- 
forded by the background turbulence, NM^ III,0 \li), are 
expressed through the similar deviations of the second 
moments, M^(k) - M^'^k) [see Eq. ©]. The su- 
perscript (0) corresponds to the background turbulence. 
First, we solve Eqs. (|A1|) - (|A3[I neglecting the sources 
ifpl^plfj with the large-scale spatial derivatives. Then 
we will take into account the terms with the large-scale 
spatial derivatives by perturbations. We assume that 
r/k 2 <C r -1 (fc) and vk 2 <C r~ 1 (fc) for the inertial range of 
turbulent flow, where v is the kinematic viscosity and r\ 
is the magnetic diffusion due to the electrical conductiv- 
ity of fluid. We also assume that the characteristic time 
of variation of the mean magnetic field B is substantially 
larger than the correlation time r(fc) for all turbulence 
scales. We split all correlation functions into symmetric 
and antisymmetric parts with respect to the wave number 



k, e.g., fij = fif+fif, where /, 



[Ai(k)+/ y (-k)]/2 

is the symmetric part and fj*' = [/y(k) — /y(— k)]/2 is 
the antisymmetric part, and similarly for other tensors. 
Thus, in a steady-state Eqs. lfAT ]) -(|M )l yield: 



1 



[(l + «fW+C(k)], 



l + 2ip 



(0), 



r ^[^ 3 (k) + (l + V)^ ) (k)], 



(A7) 

(A8) 
(A9) 



//, = ±(B.vMp + [g qj (k)(2P£\k)-5%) 



-g qi (-k)(2P$(k) - 5$)]B n , q - B n<q k n $p q , 



(A4) 



1% = |(B.V)$(f - [giq{k )S%+g jq (-k)S%]B n , q 

-B n , q k n ^ , (A5) 



?(2) 



c(2) 



% = ^(B.V)(f ij +h ij )+hi q (2P$(k)-6%)B. 



(A6) 



where Pfi\k) = 



hkj/e, 4p (k) 



<7-y(k) + 9ji{— k), and Bjj = Vj-Bi, the terms J\fij, 
Mhij and A/gy are the third-order moment terms ap- 
pearing due to the nonlinear terms, /y 9 = (l/2)dfij/dk q , 

and similarly for hij q and $L g ;. For the derivation of 
Eqs. (|A1|) - (|A3[) we use identities given in 23]. We take 
into account that in Eq. (|A3[) the terms with symmet- 
ric tensors with respect to the indexes "i" and "j" do 
not contribute to the mean electromotive force because 

&m — ^mji gij • 



where ip(k) = 2(rk-B) 2 



fij , hij and are solutions 



without the sources I\j , I-j and /?• . The correlation func- 
tions (k), hljQc) and (k) vanish if we neglect the 
large-scale spatial derivatives, i.e., they are proportional 
to the first-order spatial derivatives. 

Next, we take into account the large-scale spatial 
derivatives in Eqs. (|A1[) - (|A3[) by perturbations. Their 
effect determines the following steady-state equations for 
the second moments fij , hij and gij : 



#>(k) = ir(k.B)$^(k)+r/i , 



(A10) 



/§>(k) = -ir(k.B)^f' s) (k)+r4, (All) 
4 s) (k) = i r(k.B)(^ ) (k)-/ ll ( ; ) (k))+r^, 

(A12) 

where l>|f' s) = [$^ M) (k) + $^ M) (-k)]/2. The solution 
of Eqs. (TM0i-(EI2) yield 



$(M lS)(k) 



{^-^+ir(k.B)(7^-//, 



1 + 2t/> l y 

1 7« i 7 / J 



(A13) 
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Substituting Eq. (jAUs]) into Eqs. (|A"l0|l - l|A"l2) we obtain 
the final expressions in k-space for the tensors /y (k), 



^(k), 4 s; (k) and ' SJ (k). In particular 



~(s) 
1%j 



(M,s), 



2\(0) 



$(M, s )/ k x T(fc)l¥(fc)(u : 

+ fc nfc 4 2 ] - k mk S%) - 2( e + 2^)(k n] 8[H 



-k S {2) ) 



B 



(A14) 



The correlation functions /y (k) , h^' (k) and g^' (k) are 
of the order of ~ 0(V 2 ), i.e., they are proportional to 
the second-order spatial derivatives. Thus gtj + gtj is the 
correlation function of the cross- helicity, and similarly for 
other second moments. Now we calculate the mean elec- 
tromotive force £i(r = 0) = (l/2)ei nm J $£ M ' 
Thus, 



'(k) dk. 



l + 2i/> 



[7« n +ir(k-B)(4„-7^J]dk. 

(A15) 



We use the following model of the background 
anisotropic and inhomogeneous two-dimensional MHD 
turbulence: 



K U ,)(°)(k) = (u 2 )(°) W v (k)\(l-a v ) 



(fti 6j) (0) (k) = <b 2 ) (0) ^ m (/fc) ^ (1 - tr m ) 



r(2) 



(m) 



fc,A 



) + a" 7 % ( 5(k-B) 



(A17) 



Sjej, 5y is the Kronecker tensor, e,; 



where <5- 2) = <5y 
is the unit vector which is perpendicular to the plane 
of the two-dimensional MHD turbulence, fey = fcjfcj/fc 2 , 
(T u and er m are the degrees of anisotropy of the velocity 
and magnetic fluctuations of the background turbulence, 
and a v > a m , /% = B.Bj/B 2 , W v (k) = W m (k) = 
E{k)/2nk. The energy spectrum of the velocity and 
magnetic fluctuations is E(k) = k^ 1 {q — l)(k/ko)~ q , 
the turbulent correlation time is r(fc) = 2To(fc/fco) 1 ~ ^, , 
where 1 < q < 3 is the exponent of the energy spectrum, 
ko = l/lo, and Iq is the maximum scale of turbulent mo- 
tions, To = Iq/uq, uq is the characteristic turbulent ve- 
locity in the scale lo. The inhomogeneities of the veloc- 
ity and magnetic fluctuations of the background turbu- 
lence are characterized by — Vi(u 2 )^ /(u 2 )(°) and 



A 



(6) 



V l (b 2 )^/(b 2 )^. Note that ( Ui 6j)W(k) = 0. In 
Eqs. (|A16|) and (|A17|) we neglected small quadratic terms 
in the parameters A^' and A^ . 



After the integration in k-space we obtain £i = aijBj + 
bijkB jt k, where B^ k = V k Bj and 

a l3 = 2 Vt [(1 - OA*? - 6(1 - a m )A^} e ipn K$ , 

(A18) 



2 Vl 



[l-a v + e(l - a" 1 )] {enkK 



(i) 
pp 



=■•■ K {1) ) 



-2e(l - a m )s ink K, 



(i) 



2e infc ^. ) [l-^ 



-e(l-0] 



2r, T [3 pk {(a v +ea m )e 



kip 



Here /3 


= B/B, 


K (1) 
y 






/■27T 




~ Jo 1 


£>i(y) 


2^ , 

= — (v 




2/3 2 





-2ecr m £ mp (ex/3)„(ex/3) J ] 



:ify[y(a;)] da; , 



(A19) 



ay 



.^ = /J 1 ( 2 /)4 2) +U 2 (y)/3 y -, 



+ 1-1), D 2 {y) 
2{2-q) 

7 = r~ « c : 



2tt 



3-1 



equipartition 

field. For q = 5/3 the parameters 7 = 1 and c = 1/2, 
and for q = 3/2 the parameters 7 = 2 and c = 1/3. 

To determine the nonlinear turbulent magnetic dif- 
fusion rjij(B) we use an identity: r/ij — (ei kp bj kp + 
Sjkpbikp)/ '4. The nonlinear turbulent magnetic diffusion 
coefficient along the mean magnetic field, rj B , and the 
cross- field turbulent magnetic diffusion coefficient, T) ± , 
are given by: 



2r, 7 



e(l-0]*i(/3) 



V± = 2 VT [l-a v -e(l-a m )]^(P) 



(A20) 
(A21) 



where rj T — Iqu /2, the functions ^(fi), *i(/3) and 
their asymptotic formulas are given in Appendix B, = 
4 (B/B C q), B cq = \J (u 2 )(°) is the equipartition field and 
the param eter e = (b 2 )(°V(u 2 ) (0) . To derive Eqs. (fA^O]) 
and (|A21|) we used the following identities: 

^i£ijp/3p k ^ k Bj — , &i£ipk(3pj^J k B j Ayj4 , 

eiEi n p(3p k VkBj(ex0) n (ex0)j = A X A , 

where B = e x B. The nonlinear turbulent magnetic dif- 
fusion coefficients r\ B and for a v = a rn = are given 
in Sec. Ill [see Eqs. © and ©]. 

Now we determine the parameter e = (b 2 )(°)/(u 2 )(°) 
using budget equation for the evolution of the mean- 
square magnetic potential (a 2 ) [see Eq. ©]. To this end 



7 



we determine the energy of magnetic fluctuations (b 2 ) 
which is obtained from Eq. (|A8[) by the integration in 
k-space. The result is given by 



u 



2\(0) 



(i-O[i-0O9)] + c[i + ^ 



+</,( / 3)(l-a m )] 



(A22) 



where the function <f>{0) and its asymptotic formulas are 
given in Appendix B. Combining Eq. (|A22j) with the 
steady-state solution of Eq. (JSJ) we determine the param- 
eter 5(e) = 1 - a v - e(l - a m ): 



5(e) = 2 



Rm- 



B 2 



divF A 



1 - (J r 



-2Rmvl/(/3)_ 



<K/3) 
(A23) 



When B > B oq / VRm, Eq. (|A23|) yields the parameters 



1 - cr " 



1 - a" - 



divF^ 



4?7 T B 2 *(/3) 



(A24) 



Using Eqs. (|A20|) . (|A21j) and (|A24|) we obtain the non- 
linear turbulent magnetic diffusion in two directions: 



2ry T 



-(i-O*i08)] 



divF" 



divF" 



4r] T B 2 ^(f3) 1 



2B 2 



(A25) 



(A26) 



Note that there is a small range of the magnitudes of the 
mean magnetic field when there can be an anomalous 
behaviour of the nonlinear turbulent magnetic diffusion. 
At B — > 0.18 B e q the function *ff((3) changes sign (see 
Fig. 1). On the other hand, the function 5(e) changes 
sign for slightly larger value of the magnetic field B > 
0.18 B eq [see Eq. (|A23I) ]. Therefore, this implies that at 
B > 0.18 B cq the nonlinear turbulent magnetic diffusion 
can be anomalously large. The width of the range of the 
anomalous behaviour of the nonlinear turbulent magnetic 
diffusion is very small, SB ~ 1/Rrn. In this range the 
steady-state solution of Eq. {8J for B > 0.18 B oq does 
not exist. 

To determine the nonlinear effective drift velocity 
V cff (B) of the mean magnetic field we use an identity: 
y( eS ) — £k j iai j/2, which yields 



rcff 



-2t] T [(l - a v )A (v) - e(l - - m ) A(m) ]*i(/ 3 ) 



(A27) 

where the function v E'i(/3) and its asymptotic formulas are 
given in Appendix B. When B > B cq /VWm, Eqs. (|A"24|) 



and (|A27|) yield 

yoff = 



-2 Vt 



(l-OfAW-A^) 



div F J 



4rj T B 2 ^(f3) 



*i(/9) . (A28) 



The nonlinear effective drift velocity V cff (B) of the mean 
magnetic field for a v = a m = is given in Sec. Ill [see 
Eq. 



APPENDIX B: FUNCTIONS *(/?), *i(/3) AND <t>((3) 

In this section we present the functions ^((3), 
and <f){0) used in Sec. Ill: 



1 



= - / 

i r 1 



l + 2y(x) — 
ay. 



1 + 2y(x) 



dy 



[D 1 (y) + D 2 (y)}dx, 
Di{y)dx , 



m = ^ [ [2D 1 {y) + D 2 (y)]dx 
These functions for q = 5/3 are given by 
*(/?) 



/34 

^ [M(/3) - M(/3 Rm 1 / 4 )] + yL(/3,Rm) , 

(Bl) 



*i(/J) 
L(J3, Rm) 



1 

6/3 



2 {2 - [2 - 5/3 2 (l - 3/3 2 )] v/2/3 2 + 1} 



+ -L(/3,Rm) , 



V2/3 2 + l-— L(/3,Rm) 



(B2) 
(B3) 



= P 4 



In 



2/3 2 VRm+ 1 - 1 
2/3 2 VRm+ 1 + 1 



-In 



M(y) = 



^2(3 2 + 1-1 
v/2/3 2 + 1 + 1 
1 



5[1 - 5y 2 (l + 6y 2 )] 



+^(i-vV + i) 



(B4) 



(B5) 



Asymptotic formulas for the functions ^f{j3), v I / i(/3) and 
</>(/?) are as follows. For (3 -C Rm" 1 / 4 these functions are 
given by 



25 

1 - 9/3 2 + y/3 4 InRm 
1 - 3/3 2 + |/3 4 InRm 



1 - -0 1 InRm , 
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for Rm ^ 4 <C (3 <C 1 they are given by 



*(/?) = 
*iC9) = 



1 -9/3 2 + 50/3 4 |ln/3| 



1 -3/3 2 + 10/3 4 |ln/3| 



0(/3) = l-2/3 2 |ln/?| , 



and for (3 > 1 these functions are given by 



L7 

T 



*^ = -3^f 1 - 
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